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New constraints on Parity Symmetry from a re-analysis 
of the WMAP-7 low resolution power spectra 
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ABSTRACT 

The Parity symmetry of the Cosmic Microwave Background (CMB) pattern as seen 
by WMAP 7 is tested jointly in temperature and polarization at large angular scale. A 
Quadratic Maximum Likelihood (QML) estimator is applied to the WMAP 7 year low 
resolution maps to compute all polarized CMB angular power spectra. The analysis is 
supported by 10000 realistic Monte-Carlo realizations. We confirm the previously re- 
ported Parity anomaly for TT in the range 5£ = [2, 22] at > 99.5% C.L.. No anomalies 
have been detected in TT for a wider £ range (up to £ ma x = 40). No violations have 
been found for EE, TE and BB which we test here for the first time. The cross-spectra 
TB and EB are found to be consistent with zero. We also forecast Planck capabilities 
in probing Parity violations on low resolution maps. 

Key words: cosmic microwave background - cosmology: theory - methods: numerical 
- methods: statistical - cosmology: observations 



1 INTRODUCTION 

The anisotropy pattern of the cosmic microwave background 
(CMB), measured by the Wilkinson Microwave Anisotropy 
Probe (W MAP), probes cosmo lo gy with unprecedented pre- 
cision (see lLarson et al.l (|2010l ); iKomatsu et~aH <|2010h and 
references therein). WMAP data are largely consistent with 
the concordance A cold dark matter (ACDM) model, but 
there are some interesting devi ations from it, in particu- 
lar on the largest ang ular scales l|Copi et alj|201oh . See also 
(|Bennett et al.ll2010l ) for a critical point of view upon the 
subject. 

A large number of papers dealing with these anomalies 
have been published in the last years. We list below those 
that are the most studied: a) lack of power on large angular 
scales. The angular correlation function is found to be un- 
correlat ed (i.e., consistent wit h zero) for angles larger than 
60°. In (|Copi et al.ll2007l . 120091 ). it was shown that this event 
happens in only 0.03% of realizations of the concordance 
model, b) Hemispherical asymmetries. It is found that the 
power in temperature coming separately from the two hemi- 
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sph eres (defined by the ecliptic plane) is unlikely asymmet - 
ric (|Eriksen et al. 2004; H ansen. Bandav and Gorskill2004l ). 
It has been confirmed in the WMAP 3 year and 5 year re- 
lease (|Eriksen et al ] |2007l : lHansen" et al. 20091: iHoftuft et al.l 
2009) and it is present in the COBE data as well, al- 
though with lower significance. The temperature power spec- 
tra of the opposing hemispheres are inconsistent at 3a 
to 4cr depending on the range of multipoles considered. 
The asymmetry has been detected in low resolution maps 
|Eriksen et al.ll2004l ). both in angular and multipoles space, 
but it extends to much s maller angular scale s in the multi- 
pole range 5t = [2,600] |Hansen et all [2009). At large an- 
gular scales the hemispherical asymmetry has been tested 
for the first time in polarization maps in (|Paci et al.ll2010l ) 
making clear that this anomaly is evident only in inten- 
sity at WMAP sensitivity, c) Unlikely alignments of low 
multipoles. An unlikely (for a statistically isotropic ran- 
dom field) alignment of the quadrupole and the octupole is 
found in (Tcgmark, de Olivcira-Costa and Hamilton 2003; 
ICopi. Huterer and Starkmanl |2004|; ISchwarz et al.l |2004| ; 
Weeks 2 0041 ; lLand and Magueiidl2005al ). Both quadrupole 
and octupole ar e shown to align with the CMB dipole 
jCopi et al.ll2007tl . Other unlikely alignments are described 
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in dAbramo et all liooil; IWiaux et al. I |2006t IVielva et all 



Gruppuso and Gorskil 201oh and a test for detecting 



fore ground residuals in the al i gneme nt estimators is present 
in dGruppuso and Buriganal |2009h . d) Non-Gaussianity. 



IVielva et al. l|2004h detected a localized non-Gaussian be- 
havior in the southern hemisphere (call ed Cold Spot) us- 
ing a wavelet analysis technique (see also lCruz et al.l (2005, 
2009) ). Large scales non-Gaussian analyses can be found in 
llBernui and Reboucasll2010l;lBernui. Reboucas and Teixeiral 
l20ld ). See also lPietrobon et al.l \20M ) where the needlet for- 
malism has been applied to the WMAP 5 year data, look- 
ing for evidence of non-Gaussianity in the bispectrum of 
the needlet a mplitudes, e) Parity asym metry. It has been 
suggested in l|Land and Magueiioll2005bl ) that an estimator 
built upon the point parity symmetry might be used as a 
practical tool for detecting foregrounds. In particular these 
authors consider whether the observed low CMB quadrupole 
in temperature could more generally signal odd point-parity, 
i.e. suppression of even multipoles. However they claim that 
WMAP dataset never supports parity preference beyond 
the me agre 95% confidence level. Later, iKim and Naselskvl 
l|2010al ) found that the Parity symmetry in the temperature 
map of WMAP 3 and 5 year data is anomalous at the level 
of 4 out of 1000 in the range 81 = [2, 18]. This analysis have 
been repeated in the WMAP 7 year data confirming the 
anomaly at same level for a slightly wider range 81 = [2, 22] 
ijKim and Naselskvl[2010b1 ). 

In this paper we address the issue of the parity asym- 
metry estimating the angular power spectra (APS) of the 
WMAP 7 year data at large angular scales by an opti- 
mal and unbiased estimator such as the quadratic maxi- 
mum likelihood (QML) estimator. The same estimator is 
used to analyze 10000 simulated maps in which the noise 
is extracted from the low resolution noise covariance ma- 
trix of the WMAP 7 year data. This approach is novel since 
the estimates of the Ci are made jointly in temperature and 
polarization allowing for a global, more robust estimate of 
all six CMB spectra. In particular, it allows to extend to 
polarization an analysis that has been performed only for 
temperature so far. 

The used implementation of the QML is called BolPo l 
and it has been already adopted in jGruppuso et al.ll2009l ) 
to compute the APS of the low resolution WMAP 5 year 
data and in (|Paci et al.ll2010l ) where the hemispherical power 
asymmetries have been studied for the same data set. 

The paper is organized as follows. In Section [2] we de- 
scribe the performed analysis, introducing the property of 
symmetry we expect the CMB maps to have in Subsection 
12.11 providing the algebra of QML estimator in Subsection 
12.21 specifying the used data set and the performed simula- 
tions in Subsection 12.31 and defining the suitable estimators 
in Subsection 12.41 Results are given in Section [3] and fore- 
casts for Planck are provided in Section 0] Conclusions are 
drawn in Section [5] 



2 DESCRIPTION OF THE ANALYSIS 
2.1 Introduction 

All-sky temperature maps, T(n), are usually expanded in 
terms of Spherical Harmonics Yj m (n), with n being a direc- 



tion in the sky, namely depending on the couple of angles 
(0,4): 



a T ,im = / dD.Y e * m (h)T(h) . 



(1) 



where ar,tm are the coefficients of the Spherical Harmonics 
expansion and dfl = d8d<j> sin 8. Under reflection (or Parity) 
symmetry (n — > — n), these coefficients behave as 



-1) a T j 



(2) 



Analogously for polarizations maps, taking into ac- 
count the usual combination of Stokes parameters (Q(n) and 
U(h)) 



1±2,(m 



d£lYl xlm (h) (Q(h)±iU(h)), 



(3) 



where Y±2,im(n) are the Spherical Harmonics of spin 2 and 
d±2,lm are the corresponding coefficients, it is possible to 
show that under Parity 



where 

O-E.im = —(a-2,tm +1-2,£m)/2, 
Sfl.im = _ (d2,em — Ci-2,£m)/2i . 



(4) 
(5) 

(6) 
(7) 

Eqs. @, (®) and show that the cross-correlations 
cJ B = Cf B = 0. 

F u rther details can be found for example in l|Zaldarriagal 
Il998t) . l|Zaldarriaga and Seliaklll997h and explicit algebra is 
present in Appendix |A1 



2.2 Angular Power Spectra Estimation 

In order to ev aluate the APS we adopt the QML estimator, 
int roduced in ([T cgmark 1997) and exten ded to polarization 
in (|Tegmark and de Oliveira-Costall200ll ) . In this section we 
describe the essence of such a met hod. Further details can 
be found in (|Gruppuso et alj|2009h . 

Given a map in temperature and polarization 
x = (T,Q,U), the QML provides estimates Cf - with X 
being one of TT, EE, TE, BB, TB, EB - of the APS as: 



I' ,X' 



x'E^,x-tr(NE^,) 



where the Fjfx' i s the Fisher matrix, defined as 



dCf 



and the E x matrix is given by 

, _ i a dc ! 
Bx -2 C dd¥ c ' 



(8) 



(9) 



(10) 



with C = S(Cf) + N being the global covariance matrix 
(signal plus noise contribution). 

Although an initial assumption for a fiducial power 
spectrum is needed, the QML method provides unbi- 
ased estimates of the power spectrum contained in the map 
regardless of the initial guess, 



(C?) 



(11) 
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where the average is taken over the ensemble of realizations 
(or, in a practical test, over Monte Carlo realizations ex- 
tracted from Cf). On the other hand, the covariance matrix 
associated to the estimates, 



1\X x' 



(12) 



does depend on the initial assumption for Cf : the closer the 
guess to the true power spectrum is, the closer are the error 
bars to minimum variance. According to the Cramer-Rao 
inequality, which sets a limit to the accuracy of an estimator, 
Eq. (|12[) tells us that the QML has the smallest error bars. 
The QML is then an 'optimal' estimator. 



2.3 Data set and Simulations 

In this Section we describe the data set that we have con- 
sidered. We use the temperature ILC m ap smoothed at 9.8 
degrees and reconstructed at HealPbf] (|Gorski et al.l 120051 ) 
resolution N a id e = 16, the foreground cleaned low resolu- 
tion maps and the noise covariance matrix in (Q, U) pub- 
licly available at the LAMBDA websit e [f| for the fre q uency 
channels Ka, Q and V as considered bv lLarson et al.l (j2010r i 
for the low £ analysi s. These frequency channels have been 
co- added as follows (|jarosik et alj|2007h 



mtot = Ctot(C K a m K a + C Q 1 mq + C v 1 m v ) , (13) 

3 

C to t — C Ka + Cq + Cy 



where m%, C% are the polarization maps and covariances (for 
i=Ka, Q and V) and 



(14) 

This polarization data set has been extended to tempera- 
ture considering the ILC map. We have added to the tem- 
perature map a rando m noise realiza t ion w ith variance of 
lfj,K 2 as suggested in iDunklev et ail l|2009l ). Consistently, 
the noise covariance matrix for TT is taken to be diagonal 
with variance equal to lfiK 2 . 

We have also performed a Monte-Carlo simulations in 
order to assess the significance of our results. A set of 10000 
CMB + noise sky realizations has been generated: the signal 
extracted from the WMAP 7 years best fit model, the noise 
through a Cholesky decomposition of the noise covariance 
matrix. We have then computed the APS for each of the 
10000 simulations by means of BolPol and build two figures 
of merit as explained in the next subsection. 

Two masks are considered: one for T and one for Q 
and U. Monopole and dipole have been subtracted from the 
observe d ILC map throug h the HealPix routine remove- 
dipole (|Gorski et alj|2005r ). 



2.4 Estimators 

We define the following quantities 



a 



+/- 



(^max 1) 



+/- 

E 

i=2,e mc 



£(£+1) ax 



(15) 



1 http://healpix.jpl.nasa.gov/ 

2 http://lambda.gsfc.nasa.gov/ 



where Cf are the estimated APS obtained with BolPol for 
the power spectrum X = TT, TE, EE and BB. The sum is 
meant only over the even or odd t (and this is represented 
respectively by the symbol + or — ) with £ ma x 3. 

Therefore, two estimators can be bu ilt from Eq . (1151) as 
follow s: th e ratio R x , as performed in (|Kim and Naselskvl 
2010a) or (|Kim and Naselskvl |2uTobh . 



R~ 



C+/C* 



(16) 



and, in analo gy to what perf ormed for the hemispherical 
symmetry in (|Paci et al.ll20ld ). the difference D 



D x =Cf - C_ 



(17) 



of the two aforementioned quantities. In the following, we 
drop the index X for 7? and D specifying every time we use 
them which is the spectrum they refer to. 

For our application to WMAP data, both estimators 
have been considered for the TT spectrum but only the sec- 
ond one for the other spectra (EE, TE and BB). This is due 
the unfavorable signal-to-noise ratio of the WMAP data in 
polarization. 

For X = TB and EB we simply use the average power 



c- 



(Ir, 



E 

t=2,tmo 



1(1+1) AX 



(18) 



3 RESULTS 

APS for TT, TE and EE are given in Fig. 1 and for BB, 
TB and EB in Fig. 2. In these Figures we display both the 
BolPol estimates of the WMAP 7 year maps (black sym- 
bols) and the Monte-Carlo realizations from £ = 2 to 32 
(blue symbols). For TT and TE we display also the spectra 
provided by the WMAP team (red symbols) 0. 

In Fig. 3 we show the estimator R and D for TT aver- 
aged in 2 < £ < 22 and in 2 < £ < 33. The probability to 
obtain a smaller value than the WMAP one is 0.47% for R 
in the range 51 = [2, 22] and 3.17% in the range 8£ = [2, 33]. 
For the D estimator the probability is 0.63% in the range 
51 = [2,22] and 3.17% in the range 5£ = [2,33]. The up- 
per left pa nel of Fig. 3 recovers the same level of anomaly 
claimed in lKim and Naselskvl |2010bl ). 

In Fig. 4 we plot the percentage related to the WMAP 7 
year Parity anomaly for TT versus £ ma x in the range 10 — 40 
for the two considered estimators. As evident there is not a 
single £max for which the TT a nomaly shows up, b u t rathe r 
a characteristic scale, see also iKim and Naselskvl (|2010bh . 
For the estimator of Eq. (|16[) the percentage anomaly is 
well below 1% for almost any choices of £ m ax in the range 
[15, 25] Q. As also shown in Fig. 4 the estimator of Eq. (fT7)) 
follows closely the other estimator although it is slightly less 



3 Only these spectra are available on 
|http: //lambda.gsfc.nasa.gov/| Note that the WMAP TT 
estimates are obtained through a Maximum Likelihood code 
applied to the ILC map, whereas the WMAP TE estimates are 
obtained through a pseudo-Cf method where the foreground 
cleaned V band is considered for Temperature and foreground 
cleaned Q and V bands are taken into account for Polarization. 

4 Only for £ m ax = 21 the estimator of Eq. 1 1161) exhibits a per- 
centage which is of the order of 1%. 
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Figure 1. Estimates of TT (upper panel), EE (middle panel) and TE (lower panel) APS. Dotted lines stand for the Fiducial Power 
Spectrum, taken to be the best fit of the WMAP 7 year data. Blue lines with blue error bars represent the average and the standard 
deviation of a Monte Carlo made of 10000 sky realizations in which the full noise covariance is taken into account. Each realization 
has been analyzed with the BolPol code. Red symbols are for the WMAP 7 year estimates as provided at the following web site 
http:/ /lambda. gsfc.nasa.gov/ Note that the WMAP TT estimates are obtained through a Maximum Likelihood code applied to the ILC 
map, whereas the WMAP TE estimates are obtained through a pseudo-CV method where the foreground cleaned V band is considered for 
Temperature and foreground cleaned Q and V bands are taken into account for Polarization. Black symbols are for the BolPol WMAP 
7 year estimates where the black error bars are given through the Fisher matrix. Adopted resolution: N ail i e = 16. 




Figure 3. TT. Counts (y-axis) vs the estimator (x-axis). Upper 
histograms: Ratio for the range Si = [2, 22] (left panel) and for 
the range SI = [2,33] (right panel). Lower histograms: Difference 
for the range SI = [2, 22] (left panel) and for the range Si = [2, 33] 
(right panel) . Units for the estimator D are fiK 2 . The vertical line 
stands for the WMAP 7 year value. 
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Figure 4. TT. Percentage of the WMAP 7 year value (y-axis) vs 
/„! (x-axis). Blue line is for the ratio and the red line for the dif- 
ference. This analysis shows that there is no single £m ax for which 
the TT anomaly shows up, but rather suggests the exis tence of a 
characteristic scale, see also lKim and Naselskvl teOlObl) . 



sensitive. Therefore, we find a whole multipole range, rather 
than a single l max value, where the WMAP 7 parity anomaly 
holds. This dims significantly the case for posterior biasing. 

In Table [T] we provide the results for EE, TE and BB. 



As mentioned above, only D is considered and computed for 
the four following multipoles range 51 = [2,4], [2,8], [2,16] 
and [2, 22]. No anomalies have been found and compatibility 
with Parity symmetry is obtained. 

In Table [2] we provide the results for EB and TB where 
the estimator C is considered and computed for the same 
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Figure 2. Estimates of BB (upper panel), TB (middle panel) and EB (lower panel) angular power spectrum. Conventions as in Figure 
1. 



Table 1. Probabilities (in percentage) to obtain a smaller value 
than the WMAP 7 year one 



D 


Si = [2,4] 


SI = [2,8] 


Si = [2,16] 


Si = [2,22] 


EE 


93.09 


76.21 


44.27 


46.61 


TE 


56.35 


38.88 


24.79 


22.77 


BB 


7.97 


13.42 


11.70 


44.31 




Figure 5. EB. Counts (y-axis) vs the estimator C (x-axis). Dis- 
tribution of C for Si = [2, 22] (left panel) and Si = [2, 16] (right 
panel). Units are y-K 2 . The vertical line stands for the WMAP 7 
year data. 



Table 2. Probabilities (in percentage) to obtain a smaller value 
than the WMAP 7 year one 



c 


SI = [2,4] 


Si = [2,8] 


Si = [2,16] 


Si = [2,22] 


TB 


51.78 


39.42 


6.71 


10.55 


EB 


62.73 


69.83 


55.35 


97.70 



aforementioned four multipoles range. Both the spectra are 
well consistent with zero. Only the EB spectrum shows a 
mild anomaly in the range Si = [2, 22] at the level of 97.7%. 
This is due to five estimates from I — 18 to I = 22 that 
are systematically larger than zero. When these points are 
excluded this mild anomaly drops. For example in the range 
Si = [2, 16] the probability to obtain a smaller value than 
the WMAP one is 55.35%. The latter two estimators are 
shown in Fig. 5. 



3.1 Ka channel results 

We have considered those WMAP products that are already 
foreground reduced (see subsection I2.3|l . However in order 
to test the case of a significant foreground contamination 
we have restricted our analysis to the Ka band since such 
a channel is expected to be more polluted by synchrotron 
and free-free emission than the others. Despite the lower 
signal-to-noise ratio our analysis for the Ka channel is fully 
consistent with that of the entire dataset: no anomaly is de- 
tected in polarization, while in temperature we do confirm 
its existence at the same level. These results, in our view, 
restrict (but of course not fully exclude) the chance for a 
significant foreground contamination. The lower signal-to- 
noise ratio shows up in an increase of the standard devia- 
tions associated to the probability distribution functions of 
the estimators. For example for the Ka band data set the 
standard deviation associated to the estimator D of EE in 
the i range [2, 22] grows 3.7 times with respect to the one 
obtained with the entire data set. 
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Table 3. Standard deviation for the D estimator computed in 
the range St = [2, 22]. Units are fiK 2 . 





WMAP 7 year 


Planck 


TT 


1517.17 


1509.21 


TE 


20.19 


9.08 


EE 


0.65 


0.10 


BB 


0.69 


0.04 




Figure 6. EE. Counts (y-axis) vs the estimator R (x-axis). Dis- 
tribution of R for SI = [2, 22] (left panel) and SI = [2, 16] (right 
panel). 



the range . 



[2,22]. 


Units are fj,K 2 . 




CO 


WMAP 7 year 


Planck 


TB 


0.95 


0.19 


EB 


0.023 


0.001 



4 PLANCK FORECAST 

In this Section we take into account the white noise 
level for 143 GHz cha nnel of the Planck mission 
l|Planck Collaborationll2005l) launched into space on the 14th 
of May of 2009. As in iPaci et ail |2010l ), we consider the 
nominal sensitivity of the Planck 143 GHz channel, taken as 
representative of the results which can be obtained after the 
foreground cleaning from various frequency channels. The 
143 GHz channel has an angular resolution of 7.1' (FWHM) 
and an average sensitivity of 6 \iK (11.4 fiK) per pixel - a 
square whose side is the FWHM size of the beam - in tem- 
perature (polarization), after 2 full sky surveys. We assume 
uniform uncorrelated instrumental noise and we build the 
corresponding diagonal covariance matrix for temperature 
and polarization, from which, through Cholesky decompo- 
sition we are able to extract noise realizations. For this low 
noise level we apply the same procedure adopted for the 
Monte-Carlo simulations in Subsection 12.31 Thus, from the 
set of 10000 CMB + noise sky realizations, we find that: 



• The T based estimators (both R and D) do not change 
much since at large scale the APS for T is dominated by 
cosmic variance and not by the noise. 

• For EE, TE and BB it is possible to consider even the 
R estimator. See for example Fig. 6 where the R estimator 
is computed for EE in the range 51 = [2, 22] (left panel) and 
51 = [2,16] (right panel). 

• The standard deviations for the D and C are evaluated 
in Table|3]and|4]for 51 = [2, 22] and compared to the WMAP 
7 year ones. 



5 DISCUSSIONS AND CONCLUSIONS 

The Parity symmetry of the CMB pattern as seen by WMAP 
7 year is tested jointly in temperature and polarization at 
large angular scale. We confirm the previously reported Par- 
ity a nomaly for TT in the ran ge 51 = [2,22] at > 99.5% 
C.L. dKim and Naselskvl lioiObh . Our resolution allows to 
extend the angular range up to I = 40 (see Fig. 4) find- 
ing a decrease of the anomaly for such a wider interval of 
multipoles. No violations have been found for EE, TE and 



BB which are tested here for the first time with a method 
analogous to the one employed for TT. The cross-spectra 
TB and EB are found to be well consistent with zero. The 
analysis has been performed through a Monte-Carlo made 
of 10000 sky maps in which the CMB maps are extracted 
from the WMAP 7 year fiducial model and the noise maps 
are obtained by Cholesky decomposition of WMAP 7 year 
noise covariance matrix for Ka, Q, V channels after fore- 
ground cleaning. WMAP 7 year maps and each of Monte- 
Carlo simulations have been analyzed with BolPol that is 
our F90 implementation of the QML method. As a byprod- 
uct the full APS of WMAP 7 year data is provided at large 
angular scale, see Figs. 1 and 2. 

We have also forecasted Planck capabilities in prob- 
ing Parity violations. Considering 10000 Monte-Carlo simu- 
lations we have evaluated the improvement with respect to 
WMAP 7 year sensitivity. These are shown in Table 3 and 
4. 

It is still unknown whether this anomaly Q comes from 
fundamental physics or whether they are the residual of 
some unperfectly removed astrophysical foreground or sys- 
tematic effect. As shown in Section [4] Planck data will al- 
low one to build more precise estimators thanks to the high 
sensitivity of its instruments and the wide frequency cover- 
age. Moreover Planck is observing the sky with a totally 
different scanning strategy with respect to WMAP, which is 
a benefit from the point of view of systematic effects anal- 
ysis. Thus Planck data are awaited with great interest in 
order to confirm or discard these anomalies, making real the 
possibility to have more stringent constraints on the ACDM 
model. 
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5 Of course similar arguments apply to the other anomalies as 
well. 
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APPENDIX A: PARITY SYMMETRY 

In this appendix we explicitly show the Eqs. ([2]), ((4]) and 
©■ 



Al Temperature 

In analogy to the definition given in Eq. |T]) we write 



Replacing Eqs. (|AT"2| and (fA"T0|) in Eq. CM]) we find 



«tL = / dnY e * m (n)T^ p \n). 



(Al) 



where T^ p \h) — P[T(n)} with P being the parity operator. 
By definition P\T(n)] = T(—h), i.e. under reflection sym- 
metry h — s> — n|j that in polar coordinates is equivalent to 



9 -> 7T-0, 

(f> — > (/> + TV . 

Hence Eq. (|A1|I can be rewritten as 

a T P L = / eaiy? m (-A)T(A), 



(A2) 
(A3) 



(A4) 



where the integration variables have been changed following 
Eqs. (IA2[) . (|A3[l in order to absorb the minus in the argument 
of TQ The Spherical Harmonics Yi m are related to Legendre 
functions P™ as 



21 + 1 (£-m)\ 



P™ (cos 6*) e 



i tn d> 



(A5) 



47T (£ + m)\ 
therefore 

Y lm (-K -0,0 + *) = {-l) m {-l) e+m Y em {6, 0) , (A6) 
where it has been used that 

P f m (cos(7r - 6)) = (-l) (f+m) Pr(cos(6»)) . (A7) 
Replacing Eq. (|A"6| in Eq. (|A"H gives 

a rL = (-1)' / rfn Y L{n) T(n) = (-1)* a T , £m , (A8) 



that shows Eq. ([2]). 
A2 Polarization 

Let start considering the E mode. By definition 

a ( p 1 --IfV-P) +a (P) fA9) 

a E,tm — 2 I 2,£m ^ "-2,671 ) ' \^ ) 

where 

a ±2, em = JdQ YZ 2 ,t m (n) [Q ± iU] (P) (ft) . (A10) 

Since Q ± ill is isomorphic to a bi-dimensional vector (see 
footnote 4) than 

[Q±i(/] (P) (n) = Q (P) (n)=FiC/ (P) (w) (All) 
= Q(-fi) =f it/(-n) . (A12) 



6 This holds in 3 dimensions. Note that a Parity transformation 
in two dimensions is {x,y) — ¥ (x,—y) or (x,y) —y (—x,y). In 
other words the transformation (x, y) — > (—x, —y) is not a parity 
transformation since is equivalent to a rotation. 

7 J dQ. is invariant under Parity transformation. This is trivial 
to show in cartesian coordinates. 



dn [Xt lm {-h) Q(h) - Xle m (-n)iU(h)] (A13) 



where we have changed the integration variables considering 
Eqs. (|A2I) . (IA3|) and where we have used the conventions of 
( Zar darriagalll99a ) defining Xi/2,i m as 

A Mm (f9, <t>) = -\ <t>) + Y- a ,t m (B, <t>)) , (A14) 

X 2 , lm {6A) = ~\ {Y 2 ,t m (e,4>)-Y-2,t m (9,<t>)) ■ (A15) 



Since in l|Zaldarriagal 1 19981 ) X 1 / 2 ,im, are expressed in terms 
of Legendre polyniomials it is easy to show that 



x Mm (7r -o,<f>+w) = (-iyx lfim (e, 0) 

and 

A 2 ,, m (7r - 6, 4, + tt) = (~l) e+1 X 2 j m (6, 



(A16) 



(A17) 



Replacing Eqs. (|A16|l . (IA17|l back into Eq. (|A13[) and using 
again Eqs. (1A14|) , (|A15[) we find Eq. © 

«L P L = (-l)V^m ■ (A18) 

Repeating the same steps but starting from 



1 (j p ) J p ) 



i (P) - - —[„ 



(A19) 



(A20) 



one finds Eq. JSJ, i.e. 
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